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Abstract 

We report a new type of restoration of macroscopic isotropy (homogeniza- 

tion) in fractals with microscopic anisotropy. The phenomenon is observed in 

various physical setups, including diffusions, random walks, resistor networks, 

and Gaussian field theories. The mechanism is unique in that it is absent in 

uniform media, while universal in that it is observed in a wide class of fractals. 
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In this letter, we report a new type of restoration of macroscopic isotropy (homogeniza- 
tion) in fractals with microscopic anisotropy. The phenomenon is unique in that it is absent 
in uniform media, while universal in that it is observed in a wide class of fractals. We sus- 
pect that the phenomenon is universal enough to be observed experimentally, for example, 
in spin systems close to critical points and various transport phenomena in fractal media. 
We first discuss the Sierpihski gasket as an example of finitely ramified fractals, where the 
calculations can be performed explicitly. We then turn to the Sierpihski carpet, an infinitely 
ramified fractal, and report on rigorous results. We conclude by discussing an intuitive 
picture of the mechanism. Some results of numerical calculations are also presented. 

We note that when we discuss 'isotropy' for a deterministic regular fractal, we mean 
invariance with respect to (discrete) rotations which respect the structure of the fractal. 

Resistor network on Sierpiflski gasket. In order to illustrate the phenomenon of isotropy 
restoration, we first concentrate on the simplest example of anisotropic resistor network on 
the Sierpihski gasket, a typical finitely ramified fractal. Let n be a non- negative integer, 
and put O = (0,0), a n = (2 n ,0), and b n = (2 n-1 , 2 n ~ 1 V3). Consider the ra-th generation of 
the (pre-) Sierpihski gasket, which is a triangle AOa n b n with self-similar internal structure 
composed of triangles of side length 1, as illustrated in (Fig. [I]). Each internal vertex has 4 
bonds of unit length attached. We associate a resistor of resistance 1 with each bond parallel 
to the x-axis, and a resistor of resistance r > 1 with the remaining bonds. By repeated use 
of the star-triangle relations (Y—A transforms), this n-th level network can be reduced to 
a simple triangular network (an effective network), with resistances R£(r) in the horizontal 
bond Oa n and R^(r) in the bonds Ob n and a n b n . By definition, R${r) = 1 and R V q{t) = r. 



H n (r) measures the effective anisotropy of AOa n b n composed of resistance elements with the 
basic (microscopic) anisotropy parametrized by r = H (r). Using the star-triangle relations 
we obtain the following recursion relations for R* and R y n : 



Put 




(1) 
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RX 2i%i%(2i% + 3i%)(3i% + 2i%) 
n+1 (Rf + 6R*R% + 3Rf) (R* + 2Rl 

R\ 



n+1 R* + 2R y n " 

We see from these formula that in the anisotropic regime (H n (r) ^> 1), the effective resis- 
tances satisfy the scaling behavior 

R* n+1 (r) » 2R* n (r) , R y n+1 (r) « (3/2)i2»(r) , (2) 

while in the isotropic regime (H n (r) « 1), we have 

R x n+1 (r) » Rl +1 (r) » (5/3)^(r) . (3) 

We also see that H n (r) in ([[]) satisfies if n+ i(r) -1 = /(iJ n (r) _1 ), where 

/(x) = (4x + 6x 2 )/(3 + 6x + a; 2 ). (4) 

In particular, we see the restoration of isotropy, 

limtf n (r) = l. (5) 

n— >oo 

Fig. [| gives the calculated behaviors of the effective resistances. We see a clear signal of 
the two scaling regimes (|J) and ([|). Using @, we can calculate the rates of restoration of 
isotropy. In the anisotropic regime, we have H n+ i(r) « {3/4)H n (r), while in the isotropic 
regime, we have H n+ x(r) — 1 « (4/5)(H n (r) — 1). We can also calculate the scaling limit 
F(z) = Jim / n ((3/4)V) = z- (3/2)/ + (39/14)/ + • • •, where f n is the n-th iteration of /. 
For large r and large n(l<?i< 0(log(r)/log(4/3))) we have H~ x {r) « F((4/3) n /r). We 
can prove by standard methods using (ffl) that the scaling limit exists and that F is complex 
analytic in a neighborhood of z = 0. 

We can generalize the above consideration so that the resistors parallel to Ob Q and 
ao&o have different values. If we denote the effective resistances parallel to Oclq, Obo, 

K + R^-\- R^-\- R^) R^ {R^-vR!^ 

( K + R b n + R C n){R a +R b n+R C n) 



a b , by R a n , R b n , R c n , respectively, we find R a n+l = ^f^^X^lf&P , where K 



2(R a R b +Rb r"+r" R a ) ■ Corresponding formula for R n+1 and -R ra+1 are obtained by cyclic 



permutations of the suffixes. Restoration of isotropy lim R n jR a n = lim R c n / 'R® = 1 can be 
proved in the generalized situation. 

Restoration of isotropy is not observed in uniform media. To see this, consider a resis- 
tor network of regular square lattice, whose horizontal (resp. vertical) bonds are resistors 
of resistance 1 (resp. r). The ratio of the effective resistances for n x n size network in 
vertical direction to horizontal direction is easily seen to be r, independently of n. Thus the 
anisotropy for the resistor network of regular lattice is independent of scale. The restoration 
of isotropy which we observe on the Sierpihski gasket is a feature absent on uniform media. 

Related models on Sierpinski gasket. We described restoration of isotropy in terms of 
resistor networks The phenomenon is also observed in various other physical setups, 

including random walks and diffusions [[],f|] and Gaussian field theories || . A related math- 
ematical problem of the construction and uniqueness of diffusions on the Sierpihski gasket 
is dealt with in ||. We also remark that there is another aspect in homogenization, that 
a diffusion with microscopic irregularity restores macroscopic uniformity, as studied in |7j 
for finitely ramified fractals. This aspect, in contrast to what we deal with here, is not 
specific to fractals and has been known in Euclidean spaces. (For other related references 
in mathematics literature, see the references in ||.) 

Restoration of isotropy on Sierpinski carpet. The finite ramifiedness of the Sierpihski 
gasket implies that the recursion relations are finite dimensional, and the analysis can be 
made explicitly. One might then wonder if the isotropy restoration we found above is a 
special feature of models on finitely ramified fractals. In |9j we have proved a mathemat- 
ical theorem for a class of infinitely ramified fractals, which establishes that the isotropy 
restoration is a universal phenomenon. 

To state the result of |§, let n be a non-negative integer, and consider the pre- Sierpihski 
carpet F n , which is a subset of a unit square [0, 1] x [0, 1] obtained by removing small squares 
recursively as for constructing the Sierpihski carpet [|10] , until squares of side length 3 _n are 
reached, where we stop so that smaller scale structures are absent (Fig. Let r > 1, and 
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assume that F n is made of a material with a uniform but anisotropic electrical resistivity, 
such that for a unit square made of this material, the total resistance is 1 in the x-direction 
and r in the y-direction, and the principal axes of the resistivity tensor are parallel to the x 
and y axes. Equivalently, we assume that the energy dissipation rate per unit area for the 
potential (voltage) distribution v(x,y) is (||) 2 + \ (f^) 2 - (Note that by linear transform in 
coordinate y' = y\/r, the formula becomes that of isotropic material. Hence, in experimental 
situation, one may as well start with a rectangle made of isotropic material, with rectangular 
holes.) 

We introduce the effective resistance R^{r) °f Fn m x direction, the resistance observed 
when we apply voltage between two edges x = and x = 1. Likewise we define R^(r) and 
introduce the effective anisotropy H n (r), as in ([!]). Ho(r) = r parametrizes the anisotropy 
of the material composing F n . We can prove the following ||. 

Theorem 1 — There is a finite constant C > 1, independent of r and n, such that for any 
initial anisotropy r > 0, we have the weak restoration of isotropy (weak homogenization) 
in the sense that 1/C < H n (r) < C holds for sufficiently large n. (How large n should be 
depends on the value of r.) 

We believe that C can be taken arbitrarily close to 1, as in (El), but this is still beyond the 
reach of present mathematical techniques, for the infinitely ramified fractals. We emphasize 
that we have concrete rigorous results as Theorem 1, in spite of the difficulties for the 
infinitely ramified fractals. 

Analogous results hold if we consider a cross-wire network G n defined by replacing each 
smallest size square of F n by a horizontal and vertical cross-wire of four resistors (connected 
at the center of the square), whose resistances are 1/2 in horizontal direction and r/2 in the 
vertical direction. The results stated above for the board F n also hold for the network G n . 

Ideas for a proof of the Theorem. Theorem 1 is proved by decomposing the problem into 
the isotropic regime and the anisotropic regime. For the isotropic regime, an extension (to 
anisotropic case) of a deep renormalization group-type analysis of effective resistance for the 
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isotropic Sierpiriski carpet P,|TT| is applied, while for the anisotropic case, renormalization 



group-type picture in the neighborhood of degenerate fixed points U0,[§] holds. One of the 
key observations for the proof of Theorem 1 is that if H n (r) is very large (in the anisotropic 
regime), then H n (r) follows a scaling behavior. We can prove 

Theorem 2 — The limits lim s _1 liminf H n ((9/7) n s) = lim s _1 limsup H n ((9/7) n s) exist. 

This result says that while s = (7/9) n r and n are large, H n (r) decreases like c(7/9) n r. 
We can prove these Theorems by giving bounds controlling the n dependence of the effective 
resistances Q. Roughly speaking, we can show that in the anisotropic regime (H n (r) ^> 1), 

K +1 (r) « (3/2) i£(r) , R y n+l {r) « (7/6) J2»(r) , (6) 

while in the isotropic regime (H n (r) ~ 1), i?^ +1 (r) rs pR*(r), and i?^ +1 (r) « pR^(r). Here 
p = 1.25148 ± 1 x 1CT 5 is the growth exponent for the effective resistance in the isotropic 
case r = 1 P,lT[ . 



Based on these results, we conjecture that (H) holds also for the Sierpihski carpet, and 
that Fig. |2| schematically gives the behaviors of R%(r) and R^{r). 

Discussions. Our mathematical results are not very sharp numerically; we can only say 
that 10~ 10 < H n (r) < 10 10 , for large n. Numerical calculations for the Sierpihski carpet may 
therefore be of interest. We give results for the resistor network G n . Obviously, Rq{t) = 1 
and Ro(r) = r. It is not difficult to find Rf(r) = (3r + 4)/(2r + 3). The exact result for 
n = 2 is 

_ 324r8 + 3960r 7 + 17169r 6 + 37077r 5 + 44639r 4 + 30842r 3 + 11900r 2 + 2325r + 174 
2 ^ ~ 144r 8 + 1924r 7 + 8850r 6 + 20052r 5 + 25146r 4 + 17976r 3 + 7128r 2 + 1422r + 108 

Note that R^{r) = rR®(l/r), with which we can calculate -R^(r) and H n (r) from these for- 
mula. We have numerical results for 3 < n < 7, obtained using Gaussian relaxation method 
(Table |). We see that as r is increased, the n dependence of R%(r) rapidly approach (3/2) n , 
and that for large n, those of R^Jj) approach c(7/6) n with c = 6/5. These observations are 
consistent with (|J), implying scaling behavior in the anisotropic regime. (Deviation from 
scaling of R y n for small n in the data can be explained if we notice that we are calculating 
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the network G n instead of the board F n .) In particular, we see that for any value of r > 1, 
R^/R® monotonically decreases as n is increased, which indicates the tendency of restoration 
of isotropy. 

We expect that the scaling limit 

z lim H n ((9/7) n /z) = c + dz + --- 

n— >oo 

exists, where c is the limit in Theorem 2. The data and the fact that -R^(r) is a rational 
function of r makes it possible to find an estimate 

J5(0) = kr^jr) = c(7/6)" - 3~"/5 

with c = 6/5. Thus the constant term c in the scaling function is determined. We need 
more data to determine d, but the calculations become rapidly time consuming as n or r is 
increased. 

Let us discuss general intuitive picture of the restoration of isotropy, in terms of random 
walks The fractals may be regarded to have obstacles or holes in the space, when 

compared to uniform spaces. Intuitively, a random walker that favors horizontal motion 
performs a one- dimensional random walk between a pair of obstacles, and eventually is 
forced to move in off-horizontal direction before he could move further horizontally. There 
are obstacles of various sizes, separated by distances of the same order as their sizes, hence 
globally, the random walker is scattered almost isotropically. On uniform media such as 
regular lattices or Euclidean spaces, these obstacles are absent, hence the anisotropic walk 
keeps anisotropy asymptotically. 

The Sierpihski gasket and the Sierpihski carpet have exact self-similarity, and one may 
doubt the 'extrapolation' to figures without exact self-similarity. However, we can prove 
that the restoration of isotropy occurs for anisotropic diffusions on the scale-irregular abb- 
gaskets, a family of fractals which are scale-irregular, i.e. do not have exact self-similarity 
Qj. These considerations suggest that the restoration of isotropy is to be observed on a wide 
class of random media. For example, numerical calculations on the percolation clusters may 
provide interesting observations. 
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FIGURES 



FIG. 1. Pre-Sierpihski gasket. 

FIG. 2. R%(r) (lower plots) and R^(r) (upper plots) on the pre-Sierpihski gasket for r = 100. 
The lines are the scaling predictions @ and (||). 

FIG. 3. Pre-Sierpihski carpet F3. 
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TABLES 



TABLE I. Effective resistances -R^(r) and R%(r) for the pre-Sierpihski carpet network G. 





r = 10 


r = 100 


r = 1000 


r = 10000 


r = 100000 


n 


KM 


Rl{r) 


K(r) 


Ry(r) 


Kir) 


Kir) 


Kir) 


Kir) 


Kir) 


Kir) 


3 


2.831057 


19.64149 


3.238145 


190.6445 


3.356806 


1899.017 


3.373110 


18982.35 


3.374810 


189815.7 


4 


3.798415 


23.47825 


4.614455 


224.0274 


4.963201 


2223.085 


5.049858 


22209.29 


5.061194 


222070.3 


5 


5.070868 


28.10055 


6.524220 


263.1750 


7.258880 


2598.702 


7.524180 


25934.86 


7.585124 




6 


6.742934 


33.69136 


9.185975 


309.3891 


10.56635 


3037.488 


11.14879 


30272.61 


11.34244 




7 


8.933314 


40.46672 


12.88375 


364.0724 


15.34037 
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Fig. 3 




Fig. 2 




